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ABSTRACT 

Pulsation is believed to be the leading cause of dusty mass loss from Asymp- 
totic Giant Branch (AGB) stars. We present a temporal study of T Cep, a 
long-period Mira variable, using seven ISO SWS spectra, covering a 16-month 
period over a single pulsation cycle. The observed spectral dust features change 
over the pulsation cycle of this Mira. In general, the overall apparent changes 
in spectral features can be attributed to changes in the dust temperature, re- 
sulting from the intrinsic pulsation cycle of the central star. However, not all 
feature changes are so easily explained. Based on direct comparison with labora- 
tory spectra of several potential minerals, the dust is best explained by crystalline 
iron-rich silicates. These findings contradict the currently favored dust formation 
hypotheses. 

Subject headings: stars: AGB and post-AGB — stars: individual (T Cep) - 
(stars:) circumstellar matter — (ISM:) dust, extinction — infrared: stars 



1. Introduction 



Low and intermediate mass stars (0.8-8 M ) are major contributors to galactic chem- 
ical evolution. During the asymptotic giant branch (AGB) phase of their evolution, these 
stars lose a significant fracti on of their mass thr ough slow, massive winds at a rate of 
10~ 8 < M < 10 ~ 4 M f 7)yr~ 1 (yan Loon et al.ll2005l ). This mass loss occurs as a result of 
stell ar pulsation (IHabing|ll996l) followed by acc eleration of dust grains by radiation pressure 

Hofner et al.l 120071 ). Pulsation levitates atmospheric gas from 
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the stellar surface. Then, as the outflowing gas cools, dust can condense in the outflow 
and capture the momentum of the star's luminosity. Given that pulsation drives mass loss 
and affects the temperature and density distribution in the outflowing material, the effect of 
pulsation should be manifested in the nature of the dust, forming in the circumstellar region. 

The mineralogy (i.e. composition and lattice structure) of the dust grains is determined 
by both the prevalent chemistry, and the physical characteristics (temperature and density 
distribution) of the outflowing gas. In turn, the chemistry is determined largely by the C/O 
ratio, and to some extent by the elemental abundances in the stellar atmosphere (metallicity). 
Convection currents inside the star transport newly-formed elements (especially carbon) from 
the helium- and hydrogen-burning shells to the surface of the star. Since carbon monoxide 
(CO) is a very stable molecule and forms very easily in such stellar atmospheres, carbon 
and oxygen atoms combine into CO, leaving the more abu ndant element to dom inate the 



chemistry. Stars start their lives with cosmic C/O ~ 0.4 ( lAsplund et al.l 120091 ) and thus 



have oxygen-rich chemistry. In some cases, the dredge-up of carbon is effi cient enough to 



raise C/O ratio above unity making the star carbon-rich (e.g. IWoodl ll988). Consequently 
oxygen-rich AGB stars can have 0.4 ~C/0 < 1. 

In this paper, we investigate a single long-period Mira variable T Cep as a case study to 
determine how the mineralogy and morphology of circumstellar dust varies with the pulsation 
cycle of the central star. We present analyses of time-variations in the infrared spectra, and 
show how the spectral dust features change over the variability cycle of this Mira. In § 2 we 
discuss the history of silicate mineralogy; in § 3 we describe the observations of T Cep; in 
§ 4 we describe our analyses of the dust features, while the results are given in § 5. In § 6, 
we investigate the mineralogy of the dust and how it varies with time as well as the effect 
of grain shapes in dust spectral features. The overall discussion is in § 7, and the summary 
and conclusions are given in § 8. 



2. Previous studies: dust around oxygen-rich AGB stars 

Dust is very important in many astrophysical environments. In particular, due to the 
interplay between dust and stellar photons in driving mass loss, it is necessary to understand 
the nature of the dust around O-rich AGB stars. Dust grains absorb starlight and reemit 
photons according to their temperature and optical properties, typically in infrared (IR) 
region. Composition, lattice structure, grain size, and morphology all affect the resulting IR 
spectrum, as do the temperature and density distribution in the dust shell. We can compare 
laboratory spectra of different dust species to astronomical observations in order to identify 
the dust species present in circumstellar envelopes. 
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2.1. Finding silicates 



In the l ate 1960s, w h ile in vestigating; deviations of stellar energy distributions from 
blackbodies, iGillett et al.l (Il968h discovered a peak near 10 /zm in the spectra of four late- 
type, evolved, variable stars. Shortly thereafter a 10 /xm absorption feature w as discov- 



ered in the interstellar medium (ISM; iKnacke et al 



1969; 



Hackwell et al.l 119701 ) . Initially 



19691 ). based on s pectra of 



this feature was attributed to silicate minerals (IWoolf fc Ney 
mixtures of crysta lline silicate species predicted to form by theoretical models (jGaustad 
19631 ; lGilmanlll969l ). However, the laboratory spectra of crystalline silicate minerals s howed 



more structure within the fe ature than observed in the astronomical spectra (see IWooli 



19731 : iHuffman fc Stapplll973h. Subsequent comparison with natural glasses (obsidian and 
basal t ic glass; from e.g. IPollack et al.l Il973l ) and with artificially disordered silicates ( ID ay 
19791 ; iKratschmer &: Huffman! 1 19791 ) showed that "amorphous" silicate was a better candi- 
date for the 10 /zm feature than any crystalline silicate. 

Since the "amorphous" silicate 10 /im feature was first observed, it has been found to 



be almost ubiquitou s being found in observations of evolved stars (e.g. ISpeck et al.l 12000 



Casassus et al. 


2001 




(e.g. 


Chiar et al. 


2007) 



Casassus et al. 12001) ; m any lines of sight through the interstellar medium in our o wn galaxy 



Chiar et al.l 120071 ) : and in nearby and distant galaxies (e.g. lHao et al.l 120051 ) . Further 
spectroscopic observations of O-rich AGB stars have found that these objects show a diverse 
rang e of spectral features , which are interpreted in terms of dust cond ensation hypotheses 



(e.g. iLittle-Marenin et al.l Il990t ISloan fc Pried Il995l ; ISpeck et al.ll2000l ). In order to under- 
stand how IR dust spectra are classified and the features are interpreted we need to discuss 
the currently-favored dust formation hypotheses. 



2.2. Dust condensation sequence 



Three competing dust formation mechanisms have been considered for circums tellar 
environments: (1) thermodynamic equilibrium condensation (ILodders &: Feglevlll999l): (2) 



form ation of chaotic solids in a supersaturated gas followed by annealing (jStencel et al. 
990) ; (3) formation o f seed nuclei in a supersaturated gas, followed by mantle growth 



(jGail &: Sedlmayrl Il999l ). The latter should follow thermodynamic equilibrium as long as 



density is high enough for gas-grain reactions to occur. Several obs e rvational studies suppor t 
this thermodynamic condensation sequence (IDijkstra et al.l 120051 ; iBlommaert et al.l 120071 ) . 
however, the interpretation of the observational data is, in part, dictated by the condensation 
sequence. In 2, chaotic grains for m with the bulk co mposition of the gas, and then anneal if 
the temperature is high enough (jStencel et al.l Il990l ) . 
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Tielens 


1990; 


Grossman 


1972) 



sequence starts with the formation of refractory oxides (e.g. AI2O3) which, upon cooling, 
react with gaseous SiO, Ca and Mg to form spinel (MgAl 2 04, not shown in the figure), then 
to melilite (ranging from gehlenite [Ca2Al 2 Si07] to akermanite [Ca2MgSi20y]) and then to 
diopside (CaMgSi20e). The conversion of diopside to anorthite (CaAl 2 Si 2 08) involves a 
solid-solid reaction, and is therefore kinetically unlikely. 

In this sequence, the more abundant Mg/Fe-rich silicates form as mantles on the more 
refractory grains. Because of the relative abundance of Mg and Fe when compared to Ca 
and Al, most of the silicon atoms condense as the Mg-rich end member of the olivine 
family (forsterite: Mg 2 Si0 4 ). As the gas cools, the gaseous SiO converts forsterite into 
the Mg-rich end member of the pyroxene family (enstatite: MgSiOs). Consequently, the 
MgSi03/Mg 2 Si04 ratio increases with the decreasing temperature. At the pressure relevant 
for stellar outflows, gaseous Fe, may react with enstatite to form olivine with some iron 
([Mg, Fe]2Si04). It is also predicted that Fe condense into metallic iron grains, but it is not 
predicted that the Fe-rich end members of the olivine (fayalite: Fe 2 Si04) or pyroxene series 
(ferrosilite: FeSiOs) will form. Regardless of the precise mineral into which it condenses, 
i ron is expected to be in solid phase because it is strongly depleted from the gas phase 
JPkada et alil2009h . 



Stellar w inds are commonly assu med to be dust-driven in circumstellar environments of 
evolved stars (IFerrarotti fe Gailll2006l ). Opacity of dusty materials is an important parameter 
for the wind-driving mechanism. If iron is absent from dust, a stellar wind cannot be driven 
because of the transparency of Mg-rich silicates. Iron ha s often been include d in dust grains 
in order to increase the near-IR absorption efficiency (ITielens et al.l 119981 ). However, the 
precise effect of the iron depends on how it is incorporated, i.e. as an integral part of the 
silicate (e.g. in amorphous [Mg,Fe]2Si04 silicate) or as metallic inclusions. For Mg-rich 
silicates with metallic iron inclusions, their absorption efficiency is so high that r adiative 
heating leads to grain evaporation and only a limited amount of dust can survive (IWoitke 
2006h . 



For solar metalli city, Mg and Fe are appr oximately equally abundant in O-rich circum- 
stellar environment (ILodders fe Fegleyl Il999l ). In these dust condensation sequences, the 
condensation temperatures are dependent on the gas pressure. However, for the pressures 
relevant to the typical AGB outflow, the silicate formation temperatures are above the as- 
sumed glass transition temperature of 1000 K, implying that crystalline dust should form. 
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2.3. Spectral classification based on dust features 



Vario us attempts have been mad e to classify the observed dust features o f O-rich evolved 



stars (e.g. ILittle-Marenin et al.lll990l : ISloan fc Pried 1199a ISpeck et al.l 12000( 1 . All these au- 
thors have classified the observed spectra into groups according to the shape of the dust 
feature, which reflects a progression from a broad 9-12 /mi feature to the classic narrow 

10 /iin silicate feature. This progression of the spectral features is expected to represent the 
evolution of the dust from the early forming refractory amorphous oxides (the broad fea- 
ture) to the dominance of magnesium-rich amorphous silicates (the narrow 10 /mi feature) 
as described in § 12. 2[ 

The m o st rec ent version of IR spectral class based on silicate emission is given by 
Sloan et al.l (120031 ) . They classify the silicate dust sequence by taking the flux ratios at 10, 

11 and 12 /mi. They divided the sequence into eight segments, labeled SE1-SE8 (SE = 
silicate emission). Classes SE1-SE3 are expected to correspond to low-contrast alumina-rich 
amorphous dust seen in evolved stars losing mass at low rates and have optically thinner 
shells. Moving up the sequence, classes SE3-SE6 show structured silicate emission, with 
features at 10 and 11 /jm. The upper end of the silicate dust sequences (SE6-SE8) consist 
of sources with the classic silicate emission feature believed to be produced by amorphous 
silicate grains. These sources have optically thicker shells and higher mass-loss rates than 
sources at the other end of the sequence. 

In addition to the SE# classification, based on the features in the 9-12 fim range, an- 
other mid-IR feature at ~13/mi is sometimes present. Approximately 50% of O-rich AGB 
stars exhibit this feature and it has caused of some controversy . The so-called "13 /mi feature" 
was discovered in the IRAS LRS spectra (jVardya et al.lll986[ ) and does not lend itself to the 



standard IR feature cla ssification system. It has been found in all SE classes ( jSpeck et al. 



2000l ; ISloan et al.l 120031 ) . Furthermore, the 13 /mi f eature is associated with semiregular 



variables (SRs) ra ther than Miras or red supergiants (ISloan. LeVan fc Little-Marenin 



1996 



Speck et al.l 12000 ) . Many minerals have been proposed by various authors (e.g. ISloan et al. 
20031 : iDePew et al.l 120061 ) as the likely carrier of this feature, including corundum (crys 



talline form of AI2O3), spinel (MgA^C^) and silica (SiC^). Both corundum and spinel are 
a predicted condensate s, presolar examples of which have been found in meteorites (e.g. 
Clayton &: Nittlerl 12004 . and references therein). However, spinel has been ruled out as the 
carrier of the 13 /im feature because it also shows a strong feature at 17 /mi (irrespective of 

), and this featu r e does not occur in the spectra of AGB stars 



relative abundance of spine 
that exhi bit the 13 mm (e.g.lSloan et al.ll2003t iHeras fc Honvll2005l : IDePew et al.ll2006f ). Fur- 
thermore, DePew et al. ( 20061 ) showed that only spherical grains of spinel or corundum give 



rise to a narrow 13 /mi feature; other grain shape distributions exhibit broader features at 
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longer wavelengths. 



One of the most exciting recent findings in cosmic dust st udies was the discov ery of 
crystalline silicate dust by the Inf rared Space Observator y (ISO; iKessler et al.lll996l ) Short 
Wavelength Spectrometer (SWS; Ide Graauw et al.lll996l ). Since then, crystalline silicates 
have been found in many as trophysical environments including the solar system and extra- 
solar planetary systems (e.g. ICrovisier et al.lll996l ; lMann et al.ll2006l and refere nces therein), 



the ci rcumstellar regions of both young stellar objects (IWaelkens et al. 



1996), and evolved 



stars ( [Waters et al.lll996t iMolster et al.l l2002ai ). and in other galaxies (ISpoon et al.ll2006l ). 
These crystalline silicates were first observed around evolved stars with very hi gh mass-loss 
rates , leading to the inference that crystal formation requires such conditions (e.g. lCami et al. 



19981 ). However, recent studies on occurrence of crystalline silicates and their distribution 
with evolution suggest that cryst alline silicates preferentially occur aroun d less-evolved, low 
mass- loss rate evolved stars (e.g. ISloan et al.l 12010c iMcDonald et al.ll201lh . 



The classical amorphous silicates have broad, smooth spectral features at ~10 and 
18 /im, with no diagnostic features at longer wavelengths; that is not true for crystalline 
silicates. A cr y stalline mineralogy would manifest itself more clearly at long IR wavelengths. 
Molster et all (bo02aH c[) in their series of three papers studied dust mineralogy of several 
oxygen-rich (post-)AGB stars using ISO SWS data. By comparing the observational spectra 
with laboratory data, these authors concluded presence of Mg-rich crystalline silicates (both 
olivines and pyroxenes). However, the carrier of two of the observed far-IR features at ~20 
and 32 /im, need re-investigation. 



Various authors (e.g. iBegemann et al.l Il996l . and references therein) pointed to some 
connection between iron-rich (and possibly Mg-rich) oxide and an excess emission in the 
19-20 /xm range . Shortly th e reafte r, th e 19.5 A tm fea ture was attributed to Mg n .iFen.gO in 



depen dency by IPosch et al.l ( 120021 ) and ICamil (120031 ) . Following this study, iLebzelter et. al. 



( 120061 ) also attributed the same mineral for 20 /im feature, while studying the mid-IR dust 
features of AGB stars in the globular clusters 47 Tuc. However, the occurrence of circumstel- 
lar Mgo.1Feo.9O is not expected from thermodynamic models (e.g. lLodders fc Fegleyl Il999h 



and the low stability temperature of Mg-Fe oxides suggests that these minerals should not 
form directly from the gas phase. Consequently a large abundance of Mgo.1Feo.9O is unlikely. 
Confirmation of this attribution is further hampered by the fact that Mg-Fe oxides exhibit 
only one resonance vibrational band in the IR region, which precludes verifi cation of this 



comp osition by matching further features. A recent Spitzer Space Telescope (jWerner et al. 



2004 ) Infrared Spectrograph (IRS; iHouck et al.ll20Q4j ) study on globu lar cluster evolved star s 
also suggested that FeO is a good candidate for this ~20 /mi feature (IMcDonald et al.ll2010l ). 
In addition their modeling suggests that metallic iron grains are present around O-rich AGB 
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stars, where it is observationally manifested as a featureless mid-IR excess. 



Another interesting far-IR feature at ~32 fj,m was studied by iMolsterl (120001 ). and tenta- 
tively attributed to crystalline diopside (CaMgSi 2 6 , a pyroxene). Two recent independent 
studies of the spectrum of the AGB star RX Lac, attribu te the ~32 /im feat ure to two 



different minerals: diopside (jHony et al.ll2009l ) and fayalite ( iPitman et al.ll2010l ). Diopside 



exhibits a prominent "25 /im feature" whi ch is absent in the s pectrum of RX Lac, whereas 



all fayalite spectral features are present. iPitman et al.l (120101 ) studied the dust features of 



4 O-rich AGB stars from SE1 class (RX Lac, T Cep, T Cet, R Hya) using ISO SWS data. 
Using new laboratory data for crystalline olivines, they showed that the spectral features 
for Mg2Si04 are not present in the observed spectra. Rather the authors found that the 
peak positions and widths of the spectral features of these 4 stars match better with non- 
endmember iron-rich silicate (Mg .i8Fei. 82 SiO 4 ). The authors concluded that this result is 
unexpected, but not entirely out of the question, since it has been long a ssumed that s i licates 



in space include iron in order to account for their opacities. Moreover, iTielens et al.l (119981 ) 



found iron-bearing dust to be more refractory than typica l silicates, and therefore more likely 



to survive, which may explain the IPitman et al.l (120101 ) result. In addition to these spec- 



troscopic studies, recent work on presolar silicate grains from meteorites also suggests that 

und AGB stars, than our current models allow (e.g. 



Stroud et al. 


2008; 


Bose et al. 


2010) 



3. Observations of T Cep 



The O-rich AGB star T Ce p, was discovered by Ceraski in 1878, has been classified 

class of SE1 (i.e. 



as spectral type M5.5e - M8.8e ( lOnaka et al.lll999l ) and infrared spectra 



low-contrast broad dust feature/low mass-loss rate) with a 13 /im f eature (jSloan et al.ll2003f) . 



Its estimated distance is 210 pc from the HIPPARCOS catalogue (jPerryman et al.lll997l ). 



T Cep was observed using ISO SWS for seven times in a 16-month period (from August 
1996 to December 1997) covering one full stellar variability period. The full y processed 
post-p ipeline spectral data were acquired from an online atlas associated with ISloan et al. 



( 120031 ). Detailed data reduction information is available from the atlas webpageaj Dates 
of the observations and TDT numbers are listed in Table H], and all seven flux-calibrated 
spectra are shown in Figured] (seven observations are designated by sequential numbers). 



In addition to the IR spectral data, we have also acquired light curve data of T Cep 



http://isc.astro.cornell.edu/ sloan/library/swsatlas/atlas.html 
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from American Association of Variable Star Observers (AAVSO) databascl, which is shown 
in Figure |3l The positions in the lightcurve at which ISO SWS data was taken, are indicated 
by dashed straight lines. The apparent visual magnitudes at each ISO SWS observation are 
listed in Table [TJ According to these light curve data from AAVSO, the esti mated pulsation 



period of T Cep is 388 days, but the period shows evidence of variation. ICastelaz et al. 



(120001 ) reported the pulsation period is of 398 days, whereas, IWeigelt et al.l (120031 ) showed 
the period is of 382 days. These variation of mean pulsation period of T Cep is confirmed 
to be real by llsles fc Saw! (119891 ) and these authors also reported that T Cep vary slowly 
in amplit ude over the time. Also , an u nusual secondary variation in T Cep had been re- 
ported by Marsakova fc Andronov ( 2000 ). using the observation from AFOEV§| and VSOLJ0 



databases ([Schweitzer 



1990 



Nogamilll998l ). obtained over 75-years. These authors also re- 



ported significant variation of amplitude and asymmetry over the pulsation cycle, and they 
indicated the possibility of interference of pulsations with two periods to explain this unusual 
variation. Similar fi ndings for other several long period Mira variables are also reported by 



several authors (e.g . IWood fc Zarrolll98ll ; IZijlstra &: Bedding! 120021 ; iBenitez fc Vargas! 12002 



Ziilstra et al.ll2004 ). 



4. Analysing the dust spectral features of T Cep 

In order to investigate the potential effects of the pulsation cycle on dust formation, 
we need to consider what contributes to the observed spectra. The spectrum, Fa, can be 
interpreted as a product of the underlying continuum and an extinction efficiency factor (Q\) 
for the entire spectrum, 



F x = C x Q x x B X {T) (1) 

where B\(T) is the Planck function for a black body of temperature T, Q\ is a composite 
value including contributions from all dust grains of various sizes, shapes, crystallinities and 
compositions, and C is a scale factor that depends on the number of dust particles, their 
geometric cross section and the distance to the star. In reality the spectrum should be 
represented by: 



2 http://www. aavso.org/ 

3 ftp : / / cdsarc . u-strasbg . fr /pub/ afoe v 

4 http://www.kusastro.kyoto-u.ac.jp 
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/ a Y. x Q X}j x B^Tt) = CjXQjXB, 



(2) 



--1,3=1 



where each Bi represents a single dust (or stellar) temperature black body (of which there 
are n in total), each Qj represents the extinction efficiency for a single grain type as defined 
by its size, shape, composition and crystal structure, and each C 3 - represents the scale factor 
for a single grain type (of which there are m in total). 

For T Cep, which is a low mass-loss rate AGB star and exhibits low-contrast spectral 
features, the spectrum is dominat ed by stellar photons. Furthermore, since T Cep has an 
optically thin dust shell (see e.g. lOnaka et al.l 11999c ISloan et al.l 120031 ) . its dust spectrum 
should be dominated by the hottest and densest part of the dust shell, which is essentially 
the inner dust radius. Therefore we can simplify Eq. [5] as: 



Fx — B stellar (A, T) + BdustQdust 



(3) 



where B dust is the Planck curve for the inner dust temperature (T dust ), and Qdust is the 
emission efficiency for the innermost dust grains. This assumption is an approximation to 
the dust continuum, but the exact continuum does not significantly alter the subsequent 
analysis. 

Following Eq. [3] we can subtract the stellar contribution as approximated by an ap- 
propriate blackbody, leaving only the contribution from the dust shell. The stellar spectral 
energy distribution (SED) is assumed to be reasonably well simulated by a blackbody in the 
2500-3500 K range. The precise temperature used in our analysis, was estimated from the 
spectral type and its variability over time. For the spectral observation closest to maximum 
light (T Cep4; hereafte r T Cep,^), we u se a temperature of 3347 K, determined from the 



spectral type of M5.5 ( iPerrin et al.lll998l ); similarly at minimum light (T Cep7; hereafter 
T Cep min ) we used a temperat ure of 2566 K , deter mined from the spectral type of M8.8, 
by extrapolating the data from IPerrin et al.l ( 119981 ). For the intermediate observations we 
linearly interpolated between these minimum and maximum temperatures, relating to their 
apparent changes in visual magnitude. Table [T] lists the estimated stellar blackbody tem- 
peratures (T eff ) for all the seven observations. Figure [2] includes the original ISO spectra (in 
solid lines) together with the assumed stellar blackbodies (in dashed lines) in logarithmic 
scale (normalized at 3.0 /im, averaging over the range from 2.98-3.02 jtzmjij 



5 while there may be molecular spectral features due to water within this range, visual inspection of the 
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In reality, these spectra contain contributions of the stellar photosphere, the extended 
atmosphere, and the circumstellar dust shell. In particular, stellar atmospheric water (H 2 0) 
molecules are major absorbers in the near-and-mid-IR range. If the stellar emission is domi- 
nated by the (cooler) water layer (~2000 K), the temperatures estimated from spectral type 
are too high. However, we have repeated our analysis using stellar temperatures in range 
2000-3500 K (including ignoring the variation with pulsation phases) and found that the 
effect on spectral feature position, strength ratios and fitted dust temperature is negligible. 

Having subtracted the starlight from our spectra, we are left with the emission from 
dust (i.e. BdustQdusu Eq. |3]), which still contains a temperature factor. Therefore, we have 
fitted each starlight-subtracted spectrum with a blackbody representative of the inner dust 
temperature (T dust ). The estimated dust blackbody temperatures (T dust ) are listed in Tabled] 
and dust blackbodies are shown in Figure [2] (in dotted lines), together with the observed 
spectra and stellar continua. 

Dividing each starlight-subtracted spectrum by the best-fitted dust blackbody curve 
(normalized at 10.0 /xm, averaging over the range from 9.98-10.02 /im) leaves only the in- 
trinsic absorption/emission properties of the dust "Qdust" (shown in Figure H]), by assuming 
all dust species have the same blackbody temperature. In this case, we effectively construct 
a composite emission efficiency spectrum which can be compared directly with absorptivity 
or mass absorption coefficient measurements of minerals in the laboratory. Hereafter, we 
refer to the starlight-subtracted, dust-continuum-divided spectra as continuum-eliminated 
spectra for simplicity. 

Since, T Cep has an optically thin dust shell and its dust spectrum is dominated by the 
innermost dust, we can apply this simple modeling method, rather than the more complex 
radiative transfer modeling method. Radiative transfer modeling is known to be degenerate 
because it actually models optical depth (r), which convolves geometrical shell thickness, 
opacity and density. Furthermore, most radiative transfer models use optical constants (n Sz 
k) from laboratory experiments (or, worse, artificially derived optical constants), which are 
then applied usually using Mie theory to calculate the opacities /absorption cross sections of 



spherical grains. Re cent studies (see e.g. iMin et al.ll2003t iDePew et al.l l2006t iPitman et al. 



20081 ; ICormanl 120101 ) show that the use of spherical grains leads to unrealistic spectral fea- 
tures and this problem will be discussed further in § 16.31 In addition, mineral data in the 
form of optical constants (complex refractive index, n & k) only exist for a limited set of 



spectra shows these to be negligible, and the averaging over a range of wavelengths mitigates any potential 
problem. Furthermore since we are simply scaling a blackbody whose temperature is determined independent 
of the spectrum, small errors in scaling factor do not effect our results. 



- 11 - 



compositions. Available complex refractive indices are limited to the end-members of the 
olivine series (forsterite, fayalite), and for pyroxenes the situation is worse (only enstatite 
is available). Consequently, using radiative transfer modeling we cannot explore the min- 
eralogical parameter space so thoroughly. Although, our simple modeling approach does 
not account for different si zes, shapes and/or temperatures within the dust grain population 



sec 



Thompson et al.ll2006l . for more discussion) , it does allow us to explore the mineralogical 



parameter space comprehensively. 



5. Results from analysis of T Cep's spectrum 

Table [Tj shows that our modeled inner dust temperature (Td ust ) changes with pulsation 
cycle of the star. This result is not dependent on the precise temperature of the subtracted 
stellar black body. It is unexpected because it is usually assumed that dust forms at minimum 
light and the dust formation temperature is rather constant. However, the temperature at 
which dust can be formed depends on the mass-loss rate, which should vary through the 
stellar pulsation cycle. Thus the changing inner dust temperature may reflect the stability 
temperature at that moment in the pulsation cycle. Consequently, our simple model suggests 
that dust may not be formed constantly; rather it is sporad ic. This finding is similar to that 



for another M-type variable star Z Cyg (jOnaka et al.l 120021 ) 



We have calculated the linear correlation coefficients between the estimated stellar tem- 
perature (T e ff , which is derived from the visual magnitude and is thus a measure of the phase 
of pulsation) and the fitted dust temperature (T dust ). There is a strong correlation^ between 
them as shown in Figure EJ suggesting that the apparent change in dust temperature is 
dictated by the stellar variations, rather than the dust formation. This variation in inner 
dust temperature will be discussed further in § 16.21 

In order to study the dust properties (e.g. grain mineralogy, morphology) of T Cep, we 
concentrate on the continuum-eliminated spectra for the entire 8-45 /im, all seven of which 
are shown in Figure EJ top panel. Since, the effect of dust temperature has been removed 
in our analysis (as discussed in §H]), there is no apparent change in overall underlying slope 
over the course of the pulsation period. However, in all cases there are clear spectral features 
that peak at 9.7, 11.3, 13.1, 20 and 32 /im, indicated by dashed straight lines (see Figure |6j 
top panel). The sub-peak features within the broad 8-14 /an complex emission features (at 
9.7, 11.3, 13.1 /im) are explicitly shown in Figure [6j bottom panel. These sub-peak features 



6 The linear regression coefficient is R, but using the determination coefficient, R 2 > 0.5 is a better 
criterion for whether a correlation exists ([Thompson et al.ll2006t IChan fc Speckll201ll ) 
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are subtle, but we are confident that they are real because they do not occur at any known 
artifact wavelength of the SWS according to the SWS handbook^. 

Previous studies of the ISO SWS spectra of T Cep discovered sharp molecular bands at 
2.5, 7.3, 16.2 /im, which are attributed to molecular H 2 Q, SO-? and CO-? gas, respective ly (e.g. 
Yamamura et all Il999l ; ICami et all Il999l ; iMatsuura etaD |2002_ Ivan Malderen] l2003l l The 
attribution of the 7.3 um feature has been fu rther refined to be a combination of H2O emission 
and SiO absorption (jVerhoelst et al.ll2006l ). Figure [2] shows that these molecular emission 
features change substantially with the pulsation cycle of T Cep, whereas the dust production 
shows only slight change. Molecular emission is affected more strongly than the dust emission 
by the ste llar temperature, which changes with the pulsation cycle. The SpectraFactory 
database (ICami et al.l |2010| ) shows that molecular H 2 and OH exhibit several spectral 
bands at A > 20 /im. However, the molecular features are beyond the scope of the present 
work. 



In addition to studying the molecular features, IVan Malderen! (120031 ) also studied the 
dust features of T Cep, using the same ISO SWS spectra presented here. He concurred 
that the dust spectra of T Cep are characterized by 9.7, 11, 13 and 19.5 /xm features, but 
his interpretation of the carriers of these features is different from our present analysis; 
he used amorphous silicate (MgSiOs), amorphous alumina (Al 2 Os), spinel (MgAl 2 04) and 
amorphous magnesium-iron oxides (Mg .iFeo.gO) respectively. However, as discussed in § I2.3[ 
spinel is no longer considered the likely carrier of the 13 /mi feature, while Mgo.1Feo.9O as 
the carrier of the 19.5 /im feature also has some difficulties. Furthermore, using the same 
optical constants (see Table [2]) used by Van Malderen, it is not possible to produce sharp 
enough features to explain those observed in the spectrum of T Cep (see Figure [7]). 

In order to understand the lattice structure and the mineralogy of the dust around 
T Cep, we have sought linear correlations amongst five prominent spectral features at 9.7, 
11.3, 13.1, 20 and 32 /im. To ensure that the measure of relative strength of the features is 
independent of the choice of continuum fitting technique, we define the relative intensity as 
the ratio of the flux measured at a peak position (F pea k(Xi)) to that measured at a continuum 
point (F cont .(A 2 )). In particular the relative intensity (-F pe afc(Ai)/F con4 .(A 2 )) is less likely to be 



Tsuji et al. 


1997 


; ISpeck et al. 


2000) 



measured at Ai = 9.7, 11.3, 13.1, 20, 32 /im from each original flux-calibrated spectra. We 
used two separate continuum points (at A 2 = 8.2 and 40 /im) to verify our results. Towards 
the long wavelength end of the spectra the signal-to-noise deteriorates. The increased noise 



http://iso.esac.esa.int/manuals/HANDBOOK/sws_hb/ 
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could affect our correlations. Consequently, we take the average flux in a wavelength bin 
covering the range of 39.9-40.1 /im which contains 20 points. 

TableOlists the determination coefficients (R 2 ) for the flux ratios at peak positions (Ai = 
9.7, 11.3, 13.1, 20, 32 /tin) with respect to the continuum points at (A2 = 8.2 /im and 40 /im). 
This table shows that not all the features have strengths which are strongly/significantly 
correlated if the continuum reference is at 8.2 /mi. However, if the reference continuum point 
of 40 /im is used, almost all of the features have strongly /significantly correlated strengths. 
The discrepancy between these results for the different continuum points can be explained 
by the effect of SiO molecular absorption band, which may overlap with the 8.2 /im region 
and make this wavelength point not truly continuum. If the SiO absorption is very strong, 
then the FWHM of the feature will be affected, but not necessarily the peak positions. 

The strong correlations amongst all the peak features with respect to the 40 /im con- 
tinuum point (see Table [3]) demonstrate that the 32 /im is a real f eature (as a l so dis cussed 



in § 12.31) . rather than an artifact in these spectra, as suggested by ISloan et al.l (120031 ). The 
strong correlation between 9.7 and 11.3 /im features strongly suggests that the carrier is 
a crystalline olivine, rather than a combination of two separate miner als (i.e. amorphous 



silicate and amorphous alumina, as suggested by IVan Malderenl (120031 ) ). Furthermore the 
strong correlations between the 9-12 /im complex features and the far-IR features at (20 and 
32 /im) provide evidence that these far-IR features probably arise from the same carrier (i.e. 
crystalline olivine). Moreover, some previous studies suggest the 13 /im feature is carried 
by a different dust species from the rest of the "broad" feature (as discussed in § 12. 3p . the 
correlation of the 13 /im feature with the numero us other features supports the hypothesi s 



that this feature is due to some form of silicate (see iBegemann et al.lll997l ; ISpeck et al.ll2000l ) . 
This will be further discussed in § 16.11 The positions of the spectral features and correlations 
amongst all these features together, suggests a crystalline silicate origin. 

Having determined that the apparent broad spectral feature appears to be composed of 
overlapping features at 9.7, 11.3, 13.1, 20, 32 /im, we must now consider the implications for 
the mineralogy of the dust around T Cep. 



6. Mineralogy and morphology of dust grains around T Cep 

6.1. Comparing T Cep spectrum with laboratory data of crystalline silicates 

As discussed above, the ISO SWS spectra of T Cep, a low mass-loss rate O-rich AGB 
star, provide the evidence of crystalline silicates in its circumstellar environment, which 
call the current dust condensation sequences into question. From the strong correlations 
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among the spectral features (as discussed in § the features at 9.7, 11.3, 20 and 32/im 



Waters et al. 


1996; 


Molster et al. 


2002a: 


Pitman et al. 


2010) 



discussed in § suggest that the carrier of the 13 /im feature is a silicate, we do not have 
an exhaustive database of silicate mineral spectra to investigate this idea. In fact, silica-rich 
minerals can show a 13 /im feature and future work to find spectra of silica-rich minerals 
may be fruitful. Meanwhile, wherea s the 20 /im feature has been attributed to Mgo.1Feo.9O 



(e.g. iPosch et all l2002t ICamil 120031 ). our analysis suggests a silicate origin because of the 
correlations amongst the features. 

In order to match and identify the dust species present in the circumstellar envelop e of 
T Cep, we compare th e spect ral features of the laboratory data from lPitman et all (120101 ) and 
Hofmeister fc Pitman! (120071 ) (hereafter WashU Group) of different crystalline dust species 
to the ISO spectra of T Cep. In Figure [BJ we compare the continuum-eliminated spectra of 
T Cep max and T Cep m i n with the laboratory absorptivity data for a selection of crystalline 
olivine samples of varying composition (data taken from iPitman et all (120101 )). Likewise, 
Figure M compares the same two T Cep spectra together with laboratory absorptivity data 
for a selection of crystalline pyroxene samples of varying composition (data taken from 
Hofmeister et al. in prep.jfl 

For the olivine data, FoX is an indication of the composition such that each olivine 
has the composition Mg 2 x/iooFe 2 _ 2 (x/ioo)Si04 (dataset ranging from Fo9 to FolOO are shown 
in Figure [H]); for the pyroxene data, EnX gives the composition via Mgx/iooFei_(x/i o)Si0 3 
(dataset ranging from Enl to En99 are shown in Figure[9]). Since we have compiled laboratory 
mineral spectral data from a number of sources, these sources, along with other relevant 
sample/experimental information are listed in Table [2j The WashU laboratory spectra are 
in the form of absorbance (a), which is proportional to the optical depth, whereas we need 
to compare to the absorption efficiency (Q a bs)- 111 order to convert the laboratory data to an 
appropriate form we use absorbance (a) is proportional to the log of the absorptivity (A), 
i.e. a oc \n(A). 

In general, we compare the spectral features of the laboratory spectra of different dust 
species to the astronomical observational data in order to match and identify the dust 
species present in circumstellar envelopes. However, the spectral feature parameters (po- 



8 We have a longstanding collaboration with Anne Hofmeister at WashU and can access new labora- 
tory data through: http://galena.wustl.edu/~dustspec/info.html, even before publication. Thes e crystalline 



olivine and pyroxene serie s spectra agree well with the laboratory data from Kyoto group (e.g. iKoike et al 



2003t IChihara et al.ll2002r ) but provide fine grid spacing in composition space. Comparison among the data 



from several laboratory groups will be further discussed in 
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sitions, strengths and widths) are significantly influenced by three parameters (composition, 
temperature and grain shape of the dust grains). The effect of these three parameters are 
discussed below. 

It is clear from both Figures [8] (for the olivine family members) and M (for the pyroxene 
family members) that positions and strengths of the spectral features change with varying 
Fe/[Mg+Fe] ratio. The positions of peak features shift towards longer wavelength as the 
Fe/[Mg+Fe] ratio increases. For T Cep, the spectral features are more closely matched to Fe- 
rich silicate dust (Fo9 and Enl), rather than expected Mg-rich silicate dust (FolOO or En99; 
see Figure [8] and [9]) , which calls the conventional wisdom regarding the dust condensation 
sequence into question. These laboratory data compared to the observational spectra of 
T Cep, preclude the possibility of large abundances of the conventional Mg-rich silicates 
(FolOO, En99) to explain the dust features, unless grain shape or temperature effects can be 
invoked. 



Koike et al.l (120061 ) performed laboratory experiments to determine the effect of tem- 
perature on crystalline olivines. In general increasing temperature moves spectral features 
redwards as well as broadening and diminishing the heights of the features. This effect is 
most marked for far-IR features (at 49 and 69 /im for forsterite), and is much less effective 
in the mid-IR (8-45 /im) region. The experiments covered a wide range of sub-room temper- 
atures (8-292 K), and even this large change in temperature has little effect on the mid-IR 
features. Consequently, extrapolation of the effects to higher temperature still suggests that 
temperature is an insignificant factor in determining the positions of the spectral features 
with which we are concerned. 



Fabian et al. 


2001; 


DePew et al. 


2006; 


Sloan et al. 


2006) 



6.2. Investigation of mineralogy of dust around T Cep: Compositional 

mixtures 

It is clear from the qualitative comparisons in Figures [8] and [9] that Fe-rich crystalline 
silicates (Fo9, Enl) are promising constituents for the dust around T Cep. 

Figure [TUl shows comparison of continuum-eliminated spectra of T Cep max and T Cep m i n 
with mixture of these potential Fe-rich crystalline silicates minerals in different ratios. To 
account for the 13 /im feature, we have also included spherical corundum (A1 2 3 ) grains. We 
use this mineral simply to give a feature, and not because we necessarily believe the carrier 
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to be corundum. The source of the laboratory data of corundum (Cor) and its chemical 
composition are listed in Table [2j 

The best fit models from our calculations are mixtures of crystalline fayalite (Fo9; not- 
quite-endmember of olivine family), ferrosilite (Enl) and corundum (Cor) in varying pro- 
portion (with the error-bars of 10% for each constituent). Although the models are not a 
perfect fit, there appears to be a variation from T Cep max to T Cep m i n . This variation is best 
explained by changing the ratio of Fo9 to Enl. This is demonstrated in Figure [TDJ which 
shows the best fit models with T Cep max and T Cep min along with the laboratory data of 
each constituent separately. 

Because laboratory optical properties of solids are produced for a variety of applications, 
the various laboratory data are published in several different types of units (e.g. absorbance, 
mass absorption coefficient, complex refractive index etc.) Using the appropriate conversion 
factors, all are converted to absorptivity (A) before attempting to make mixtures for spectral 
fitting. 



For Fo9 and Enl the origin al spectra were in absorbance (a) units (e.g. Hofmeister et 



al. 



m prep; 



Pitman et al.ll2010f). which w ere converted to absorptivity (A), using A oc e° 



And for Cor, we used the iMin et al.l ( 120031 ) method to calculate the absorption cros s section 
(C a hs ) for spherical grains of size 0.1 /im using its optical constants (data taken from iGervais 
199ll ). Absorptivity (A) is directly proportional to the absorption efficiency Q a bs of a grain, 
which, in turn, is directly proportional to the absorption cross section (C a b s )- 



We use the following linear conversion equation to convert C a b s to A for corundum: 



A = C abs xnx d (4) 

where n is the number of particle per unit volume of the material used, and d is the path 
length. For thin film, d is taken as 1 /im. And to calculate n, we use 



M mo i x m H2 

where p is the density of corundum (4.02 gm/cm 3 ), m# 2 is 1.672xl0 -24 gm. And M mo i 
for A1 2 3 is (27x2+16x3=102). Using these values in EqEl we get n = 2.35xl0 22 cm" 3 . 
Knowing C a bs, n and d, we calculated A for AI2O3 by using EqlU 

Figure [10] demonstrates that these mixtures produce a reasonably good match to the 
detailed shape and as well as overall shape of the spectrum for both T Cep max and T Cep min . 
Consequently, we can conclude that the compositions of the individual dust constituents 
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remain the same in both spectra, while the relative amounts of these constituents may 
change, but are consistent with no variations at all. The spectra are consistent with a 
variation in the ratio of Fo9/Enl such that it is doubled at maximum light (T Cep max ) and 
the ratio is halved at minimum light (T Cep m i n ), while the relative abundance of Cor remains 
unchanged. This suggests that the olivine grains are slightly favored at maximum light. 

The apparent subtle changes in mineralogy with pulsation cycle may be entirely due to 
statistical effects. However, in the classic condensation sequence (shown in Figured]), Mg- 
rich olivine exists at higher temperatures than Mg-rich pyroxene. Furthermore the Mg-rich 
olivine reacts with SiO gas to form Mg-rich pyroxene. If a similar reactive process exists 
for the iron-rich endmembers, the occurrence of higher olivine/pyroxene compositions at 
maximum light, may be a temperature effect. In this case, the Fe-rich pyroxene (Enl) is 
more easily destroyed than the Fe-rich olivine (Fo9) at maximum light; while at minimum 
light the reaction with SiO gas is promoted. Essentially, the change in stellar radiation field 
causes selective destruction/processing of the inner most grains. 

From our analysis, the most striking result is that we need Fe-rich crystalline sili- 
cates in order to explain the spectral features. This result is unexpected, but supported 
by recent studies of silicate presolar grains from meteorites which show that silicate grains 



non-e quilibrium formation processes (e.g. 



Stroud et al. 


2008; 


Vollmer et al. 


2008; 


Bose et al. 



2010h . 



6.3. Investigation of morphology of dust around T Cep: Grain shape effects 

Our analysis on dust compositions around T Cep strongly suggest that the crystalline 
silicates are almost completely Fe-rich, with little evidence of Mg-rich silicates. This result 
conflicts with both current dust formation hypotheses, and studies of cosmic crystalline 
silicates to date. However, as mentioned in § 16.11 it is possible that grain shape effects may 
allow Mg-rich silicates to match the positions and strengths of the observed dust features of 



on their grain shapes (see 


Bohren & Huffman 


1983; 


Bohren et al. 


1983 


Min et al. 


2003; 


Sloan et al. 


2006; 


DePew et al 


2006 


IrCoike et al.l < 


>010 


) 



the grain shape effects of crystalline silicates. 

The spectral features of crystalline silicate grains exhibit wide variety in the positions, 
widths and strengths of their peaks. For a given series, these are dependent on the com- 
position (Fe/[Mg+Fe] ratio) of the silicates, as well as the temperatures (as discussed in 
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§ 16. ip and shapes of dust grains. Disentangling these competing effects is difficult because 
we do not have the data to test both grain shape and compositional effects simultane- 
ously. There are several laboratories that have produced spectral data for crystalline sili- 



cates in the olivine and pyrox ene series. The three 
these studies are WashU 



Chihara et al. 



(e.g. 



2002 



Fabian et al. 



Hofmeister fc Pitman 



Koike et al. 



2001 



2003, 



2006 



Jager et al.lll998l ) groups. 



aboratories, which are m ost noted for 
20071: IPitman et alJboioh. K yoto (e.g. 



Murata et al. 



2009 



Koike et al. 



20101 ) and Jena 



A series of opacity (mass absorption coefficients) measurements for the olivine series, 
covering several Fe/[Mg+ Fe] ratios from for sterite (FolOO) to fayalite (FoO) has been pub- 
lished by the Kyoto group (IKoike et al.ll20031). A finer grid of compositions across the olivin e 



series was published by the WashU group ( IHofmeister &: Pitman 



2007 



Pitman et al.ll2010h 



which agree with the Kyoto data. Here we use to finer grid of IPitman et al.l (2010 ). A similar 
series of measurements for pyrox ene from enstatite (E nlOO) to ferrosilite (Enl) has been also 
published by the Kyoto group (IChihara et al.ll2002l ). Here, we use currently unpublished 
pyroxene data from the WashU group (Hofmeister et al. in prep.), which also agree with the 
Kyoto data. 

While the data from WashU and Kyoto groups provide a good sampling over several 
Fe/[Mg+Fe] ratios for both olivines and pyroxenes, they do not measure the complex indices 
of refraction or dielectric constants for all these compositions. As a consequence we cannot 
use them to analyze the effect of grain shape on dust spectra. Thus we cannot fit both 
the compositional effect (Fe/[Mg+Fe] ratio) and the shape distribution simultaneously. For 
this reason, extensive studies of grain shape effects have generally been limited to a single 
composition: forsterite (FolOO), supplied by the Jena group. 

The complex refractive indices for the end-member compositions of olivine (Mg-rich 



Fabian et al. 


2001; 


Jaeer et al. 


1998) 



Jena grourj^ (e.g. 

for the vibrational dir ections parallel to the three crystallographic axes x, y, z. Meanwhile, 
Mukai fc Koike! (Il990f ) provided optical constants for forsterite (FolOO), but this data was 



for unoriented samples and therefore represents an average of the three axial directions. A 
list of all the sources of these laboratory data used, can be found in Table |5J 

Since solid particles in astrophysical environments are expected to be very irregular 
in shape, it is very difficult to characterize the shape of the particles in a simple way. 



9 Tcchnically Fo90 is still called forsterite and most natural samples are closer to Fo90 than FolOO 
10 optical constants available at http://www. astro. uni-jena-de/Laboratory/Database/databases. html 
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Traditionally, cosmic dust grains have been assumed to be spherical for simplicity, however, it 
is becoming incre asingly oby i ous that the use of spheri c al grains leads to un r ealistic spectra l 



features (see e.g. iMin et all 120031 : IPePew et all 120061 : IPitman et al.l 120081 : ICormanl BOlOf l. 



There are several possible approaches to addressing the grain shape effect, which we will 
compare here. 



Fabian et al.l (120011 ) used their optical constants of crystalline Mg-rich olivine (Fo90) 



and fayalite (FoO) to calculate mass absorption coefficients for IR region for a variety of 
grain shapes. These authors primarily chose: spherical grains, ellipsoids elongated along the 
z-axis, and distributions of ellipsoids; based on the assumption that mineral forms crystalline 
ellipsoids elongated along the z-axis, when condensing from the gas phase at high temper- 
ature. Whether there should be a preferred growth axis for a grain shape or not, will be 
discussed below. These authors compared their calculated mass absorption coefficients data 
directly with ISO SWS observational data of evolved stars. They concluded that the spec- 
tral features are shifted in wavelength by grain shape effects and this effect is prominent in 
mid-IR region, whereas, the spectral features in far-IR region remain practically unaffected 
by the particle shape. 



Recent experimental studies by lTakigawa et al.l ( 120091 ) on the evaporation of single crys- 
tals of forsterite (FolOO) show that this mineral evaporates anisotropically, which may lead 
to distinct grain shape distributions (e.g. disk-shaped dust grain, flattened along the y or z 
axes). They showed that forsterite always evaporates anisotropically, but the details of the 
anisotropy (i.e. which axis is elongated) depends on the experimental conditions (total gas 
pressure and temperature). They conclude that the precise shape of forsterite grains depends 
on the formation conditions in space. Because the peak positions and relative strengths of 
dust spectral features are expected to depend on grain shape, these spectral parameters may 
be diagnostic of the formation and heating conditions for the dust grains. 

Another recent study by the Kyoto grou p investigated how the spectral features of 
forsterite (FolOO) are affected by grain shape (IKoike et al.ll2010l ). They measured the IR 
mass absorption coefficients of forsterite grains of a variety of shapes including irregular, 
plate-like with no sharp edges, elliptical, cauliflower and spherical. They concluded that 
spectral features of forsterite in mid-IR region (at 11, 19, 23, 33 /xm) are extremely sensitive 
to particle shape, whereas, the features i n far-IR region (at 49, 69 jttm) remained unchanged 
despite of the different grain shapes (c.f. iFabian et al.ll200ll ). 



Min et al.l ( 120031 ) adopted statistical approach (following iBohren fc Huffman! Il983l ) in 



order to understand the effect of grain shape on IR dust features of forsterite (FolOO). They 
considered the limiting case, where the particles are very small compared to the wavelength 
of radiation (Rayleigh domain). In this statistical approach, the scattering and absorption 
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properties of irregularly shaped particles can be simulated by the average properties of a 
distribution of simple shapes (such as ellipsoids, spheroids, and hollow spheres). They cal- 
culated the absorption and scattering cross sections of different grain shapes and concluded 
that there is a strong similarity between the absorption spectra of distributions of various 
non-spherical homogeneous particles (e.g. ellipsoid, spheroids) and a distribution of hol- 
low spheres in Rayleigh domain, but that homogenous spherical particles show a markedly 
different spectrum. 



Using the approach of IMin et al.l (120031 ) and assuming that there is no tendency for a 
certain axis to be elongated, we have calculated the absorption cross-sections (C a b s ) for three 
compositions of olivine (FolOO, Fo90 and FoO; see Table [2] for sources of complex refractive 
indices) . For each composition, we calculate C a b s for four grain shape distributions: Spher- 
ical particle (SPH), Continuous distribution of ellipsoids (CDE), Continuous distribution of 
spheroids (CDS), and Distribution of hollow spheres (DHS0 The grain size is defined by the 
grain volume being equivalent to a sphere of radius 0.1 /xm. When data is available for the 
individual crystallographic axes, we calculate the C a b s separately and then average the three 
axial absorption cross sections ((C a b s ))- The final (C a b s ) spectra are shown in Figures [TTHT51 

For all three compositions (FolOO, Fo90 and FoO), the spectral feature parameters (e.g. 
peak positions and strengths) are similar for both ellipsoids (CDE) and spheroids (CDS). 
For the hollow spheres (DHS) the peak positions are similar to both CDE and CDS, but 
the spectral features are stronger. The spectral features of spherical (SPH) particles are 
significantly different in both position and strength from those of CDE, CDS and DHS. 

Figure [15] compares the positions of the spectral features in the derived absorption 
cross-sections with those observed for T Cep. It is clear that fayalite (FoO) provides a 
better match with the observed peak positions (at 20, 3 2 /xm) of T Cep, than does forsterite 



(FolOO). Following the work of iTakigawa et al.l ( 120091 ). it is possible that a single growth 



axis is elongated. However, Figure [15] shows that even using only a single crystallographic 
axis does not produce an Mg-rich olivine with features that match T Cep. 



7. Discussion 

Our analysis of the spectra of T Cep strongly suggest that the dust forming around this 
Mira is both highly crystalline and iron-rich. This has implications for not only for other 
stars exhibiting these crystalline features, but also for conventionally Mg-rich condensation 



11 Hollow spheres are meant to simulate fluffy grains 



-21 - 



sequence observed for amorphous circumstellar silicates. 



Woitkd ( 120061 ) showed that while carbon dust around an AGB star could drive a 



radiation-pressured wind, oxygen-rich dust (silicate) was too transparent. However, this re- 
sult assumes the silicates are Mg-rich. Iron-rich silicate grains tend to have higher optical/near- 
IR opacities which facilitate the capture of momentum from the star through radiation pres- 
sure. The inclusion of iron-rich silicate grains may solve this problem. 



Woitkd ( 120061 ) also showed that the dynamics in the dust-forming zones around carbon- 



rich AGB stars lead to inhomogeneous dust formation, producing fine scale structure in 
the density of the dust envelope. In these models the only condensate considered is amor- 
phous carbon. In an oxygen-rich environment, there are many potential minerals that can 
be formed and their stability is sensitive to the precise co nditions. In ad dition to the tur- 
bulent /hydrodynamic density inhomogeneities predicted by Woitkel (20061), puls ation shocks 
are expected to have a strong effects on local conditions (e.g. ICherchneffl 120061 ) . This com- 
bination of physical effects should lead to non-equilibrium dust formation and may lead to 
unexpected dust-forming conditions. Therefore we suggest that even at low mass-loss rates 
the density structure in the outflows of AGB stars is such that crystalline silicates may form, 
though it is not clear why iron-rich silicates are apparently favored. 



8. Summary & Conclusion 

We have presented an analysis of the time variations of the IR dust spectrum of optically 
thin O-rich AGB star, T Cep. 

We found that: 

1. The inner dust temperature of T Cep is variable. 

2. Dust formation is likely to be sporadic, not continuous, and has approximately the 
same composition all the time. 

3. While the observations are consistent with a constant olivine-to-pyroxene ratio, they 
can accommodate small variations with stellar pulsation. This variation in composition 
is very subtle and can be explained by selective destruction/processing of the inner most 
grains with the change in stellar radiation field. 

4. The strong correlations between the observed spectral features suggest that they all 
have same crystalline mineral as a carrier. 
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5. The structure within the broad 8-14 /xm feature with overlapping sub-features at 9.7, 
11.3, 13.1 /iin is explained by mixtures of crystalline silicates. This confirms the pres- 
ence of crystalline minerals around low-mass-loss rate O-rich AGB stars. 

6. The peak wavelength of the features at 20 and 32 fxm suggest the presence of Fe-rich, 
rather than the expected Mg-rich silicates. This can be explained as occurrence of 
non-equilibrium condensation mechanism in the outflow of the central star. 

The analysis presented here shows that our understanding of the formation of crystalline 
silicates and the inclusion of iron in those silicates is in its infancy and needs to be revised 
according to the present findings. 
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Fig. 1. — Predicted condensation sequence for O-rich environments from iTielend (119901 ) 



(after iGrossmanl Il972l ) . 



-30 - 



2 4 8 16 32 2 4 8 16 32 




1 1 — | — i — i — i | i i — i — i | i i — f 

2 4 8 16 32 

Wavelength in jam 



Fig. 2. — JiS'O SWS spectra of T Cep together with both stellar and dust continua. solid 
lines are the observed spectra; dashed lines are fitted stellar blackbody continua; dotted lines 
are dust blackbody continua (see text for details). The estimated temperatures of the stellar 
blackbodies (T eff ) and the fitted temperatures for the dust blackbodies (T dust ) are listed in 
Table ffl 
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Fig. 3. — Light curve of T Cep over a 16 month period. The seven ISO SWS observations 
are indicated by dashed straight lines. 
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Fig. 4. — Starlight subtracted, dust-continuum-divided spectra (continuum-eliminated spec- 
tra). 
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Fig. 5. — Correlation between assumed stellar temperature (T eff ) and fitted inner dust tem- 
perature (Tdust), along with the linear correlation coefficient. 
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Wavelength in |im 

Fig. 6. — Continuum-eliminated spectra of T Cep in mid-IR region. The position of the 
spectral features (at 9.7, 11.3, 13.1, 20, 32 /im) are indicated by dashed straight lines in the 
top panel. The sub-peak features within the broad 8-14 /im emission (at 9.7, 11.3, 13.1 jum) 
are explicitly shown in the bottom panel. 
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Fig. 7. — Comparison of continuum-eliminated spectra of T Cep4 (T Cep max ) with several 
amorphous minerals suggested by Van Malderen (2003). The dotted straight lines indicate 
the positions of dust spectral features as shown in Figure [6j 
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Fig. 8. — Comparison of continuum-eliminated spectra at maximum (T Cep max ) and min- 
imum (T Cep min ) light with the laboratory absorptivity data for the olivine solid solution 
with varying Fe/(Mg+Fe) ratios. FoX (X=9, 31, 50, 67, 80, 91, 100) defines the composition 
such that each olivine has the composition Mg 2 x/iooF e 2-2(x/ioo)Si0 4 . The dotted straight 
lines indicate the positions of dust spectral features as shown in Figure |6l 
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Fig. 9. — Comparison of continuum-eliminated spectra at maximum (T Cep max ) and min- 
imum (T Cep min ) light with the laboratory absorptivity data for a pyroxene solid solution 
with varying Fe/(Mg+Fe) ratios. EnX (X=l, 12, 40, 55, 99) defines the composition such 
that that each pyroxene has the composition Mgx/iooF e i-(x/ioo)Si0 3 . The dotted straight 
lines indicate the positions of dust spectra features as shown in Figure [6j 
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Fig. 10. — Comparisons of continuum-eliminated spectra at maximum (T Cep max ) and min- 
imum (T Cep m in) light (solid lines) along with the best fit model (dashed lines) of mixtures 
of different potential crystalline minerals with different ratios. The laboratory data of in- 
dividual crystalline minerals (Fo9, Enl, Cor) are also included. The dotted straight lines 
indicate the positions of dust spectra features as shown in Figure [6j 
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Fig. 11. — Absorption cross section of very Mg-rich Olivine (Mg L9 Feo.iSi0 4 ) for spheroid 
(CDS: left panel) and ellipsoid (CDE: right panel). The upper three plots correspond to 
calculated C^ s parallel to the x, y and z axes. The bottom plots correspond to average 
((Cabs))- Vertical dashed lines show positions of observed spectral features of T Cep at long 
wavelength (20 and 32 jum). 
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Fig. 12. — Absorption cross section of very Mg-rich Olivine (Mg L9 Feo.iSi0 4 ) for spherical 
(SPH: left panel) and hollow sphere (DHS: right panel). The upper three plots correspond 
to calculated C a b s parallel to the x, y and z axes. The bottom plots correspond to average 
((Cabs))- Vertical dashed lines show positions of observed spectral features of T Cep at long 
wavelength (20 and 32 /im). 
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Fig. 13. — Plot of Absorption cross section of fayalite (Fe 2 Si0 4 ) for spheroid (CDS: left 
panel) and ellipsoid (CDE: right panel). The upper three plots correspond to calculated 
C a bs parallel to the x, y and z axes. The bottom plots correspond to average ((C a b s )) of 
them. Vertical dashed lines show positions of observed spectral features of T Cep at long 
wavelength (20 and 32 /im). 
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Fig. 14. — Plot of Absorption cross section of fayalite (Fe 2 Si0 4 ) for spherical (SPH: left panel) 
and hollow sphere (DHS: right panel). The upper three plots correspond to calculated C a b s 
parallel to the x, y and z axes. The bottom plots correspond to average ((C a bs))- Vertical 
dashed lines show positions of observed spectral features of T Cep at long wavelength (20 
and 32 /im). 
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Fig. 15. — Calculated absorption cross section (C a b s ) for four shape distributions (SPH, 
CDE, CDS, DHS) of 0.1 /^m-sized grains of three different members of olivine family along 
with the directly-measured laboratory absorption spectrum. Right panel: forsterite (FolOO), 
center panel: very Mg-rich olivine (Fo91: closest compositionally in the lab spectra), left 
panel: fayalite (FoO). Vertical dashed lines show positions of observed spectral features of 
T Cep at long wavelength (20 and 32 jtim). The solid lines are the calculated absorption 
cross-sectio ns while dashed lin es (at the top) show the laboratory absorption spectra (data 
taken from lPitman et allhoioh . 
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Table 1: Information regarding observational data and blackbody curves. 



Observation 


TDT 


Observation 


Visual 


^eff 


Tdust 


Number 


Number 


Date 


magnitude 


(K) 


(K) 


T Cepl 


26300141 


08/05/1996 


9.7 


2703 


575 


T Cep2 


34601646 


10/27/1996 


9.8 


2683 


550 


T Cep3 


42602251 


01/15/1997 


8.0 


3035 


600 


T Cep4 


51401256 


04/13/1997 


6.4 


3347 


650 


T Cep5 


57501031 


06/16/1997 


7.7 


3093 


600 


T Cep6 


66101436 


09/07/1997 


10.0 


2644 


550 


T Cep7 


74602101 


12/01/1997 


10.4 


2566 


500 
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Table 2: Sources of laboratory mineral data. 



Sample 


Chemical 


Designated 


Original form 


References 




Composition 


by 


of lab data 




Amorp. enstatite 


MgSi0 3 




optical constants (n, k) 


Jager et al. (1994) 


Amor p. alumina 


A1 2 3 




optical constants (n, k) 


Begemann et al. (1997) 


Cryst. spinel 


MgAl 2 4 




optical constants (n, k) 


Fabian et al. (2001) 


Amorp. Mg-Fe oxides 


Mgo.1Feo.9O 




optical constants (n, k) 


Henning et al. (1995) 


Cryst. Olivine 


(Mg,Fe) 2 Si0 4 


FoX 


absorbance (a) 


Pitman et al. (2010) 


Cryst. Pyroxene 


(Mg,Fe)Si0 3 


EnX 


absorbance (a) 


Hofmeister et al. in prep. 


Cryst. alumina 


AI2O3 


Cor 


optical constants (n, k) 


Gervais (1991) 


Cryst. Forsterite 


Mg 2 Si0 4 


FolOO 


optical constants (n, k) 


Mukai & Koike (1990) 


Cryst. Mg-rich Olivine Mgi. 9 Feo.iSi0 4 


Fo90 


optical constants (n, k) 


Fabian et al. (2001) 


Cryst. Fayalite 


Fe 2 Si0 4 


FoO 


optical constants (n, k) 


Fabian et al. (2001) 
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Table 3: Correlations between flux ratios for th e observed spectral features. 
Continuum measured at — > 8.2/xm 40/im 
Peak Position R 2 



F9.7 


F11.3 


0.913 


0.958 


F9.7 


F13.1 


0.489 


0.899 


F9.7 


F20.0 


0.576 


0.863 


F9.7 


F32.0 


0.439 


0.579 


F11.3 


F13.1 


0.334 


0.797 


F11.3 


F20.0 


0.633 


0.826 


F11.3 


F32.0 


0.504 


0.537 


F13.1 


F20.0 


0.545 


0.869 


F13.1 


F32.0 


0.567 


0.530 


F20.0 


F32.0 


0.215 


0.323 



Coefficients of determination (R 2 ) > 0.5 constitutes a correlation; 
bold designates to which there are strong/significant correlation. 



